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Abstract 



. We study the properties of H-dibaryon matter through the SO (3) Skyrmion solution on 

a three-dimensional hypersphere S^. As the density increases, the swelling of H-dibaryon 
^ ' is found. Above a critical density, the system becomes uniform in terms of the baryon 

CO . density. In this uniform phase, the critical order parameter is largely reduced, which 

I can be interpreted as the chiral symmetry restoration. From the comparison with the 

00 I SO(2)xSO(2) solution, the S0(3) soliton is found to be the true ground state for high 

■ density system with R < 0.64 fm in the B = 2 sector with A^f = 3. 
■ 

1 Introduction 

Recently exotic hadron states have been studied elaborately. Among them, an H-dibaryon 

■ has been an object of many worksfll] as a hopeful candidate of state since the prediction by 
>• . JaffepI in the MIT bag model. Lots of theoretical studies have been made on the problem of the 

^ ! H-dibaryon, and have predicted the H-dibaryon mass below the AA threshold [|l], ^. The Skyrme 
^ \ model was first applied to the H-dibaryon by Balachandran et al. proposing the B = 2 SO (3) 
- - ' soliton configuration, which has the symmetric configuration in terms of the flavor SU(3)[|, |^. 

They obtained the H-dibaryon mass lighter than two A's (Mh ~ 1.92Mb) in the SU(3) flavor 

symmetric limit. 

Experimentally, several candidate events have been reported]^, but any conclusive result 
for its existence is not obtained yet. The binding energy of two A's in a double hypernucleus 
places a restriction on the H-mass. A double hypernucleus event found at KEK^] implies that 
the lower limit of the H-dibaryon mass is a few tens of MeV below the AA threshold. 

The possibility of the existence of H-dibaryon matter in a high-density neutron star is 
pointed out by Tamagaki et al.^. In their argument, the compactness of the H-dibaryon 
radius is an important assumption on the existence of H-dibaryon matter, although the size 
of H-dibaryon may vary in the dense system as the nucleon swelling phenomenon in nuclear 
matter [1^. If so, nuclear matter may decay into H-matter at a critical density Pc = 6 ~ 9po 



(po — 0.17 fm ) before the phase transition into quark-gluon plasma. 
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The properties of nuclear matter at high densities have been investigated in many ways. 
Castillejo et al. studied Skyrmion matter by placing SU(2) hedgehog Skyrmions on a lattice 
and show that the deconfinement phase transition and the chiral symmetry restoration occur 
above a critical density [|1T|. Manton and Ruback obtained the similar results by a simple and 
mathematically sophisticated approach. Instead of placing many Skyrmions in the flat space 
B? , Skyrmion matter is approximated by one Skyrmion on an hypersphere with radius R. 
They showed that beyond a critical density the lowest energy solution is given by the "identity 
map" and the baryon density distributes uniformly on the hypersphere [0]. In this phase, the 
spacially averaged value of the chiral order parameter becomes zero, which may imply the chiral 
symmetry restoration at high densities. This phase transition was investigated also in the case 
of two Skyrmions on S'^]!^. The result is qualitatively the same as one Skyrmion case on S"^, 
and they concluded that these behaviors are a general tendency in the Skyrme model. 

In this paper, we investigate the feature of H-dibaryon matter by using the SO (3) Skyrmion 
with baryon number 2 located on hypersphere. We study how the QCD phase transition 
occurs in this situation. We study also how the size of the H-dibaryon varies in high-density 
H-matter. For simplicity on the mathematical and physical insight, we deal with the chiral 
limit in this paper. 



2 Formalism 

The Skyrme Lagrangian in the chiral limit |T3] is given by 

C = -S-Trid^U^d^^U) + ^Tr[d^UU\ d^UU^f + £wz, (1) 

where f^^ and e denote the pion decay constant and the Skyrme parameter, respectively. Here, 
U{x) = Qi^^"-{x)/f^ ^ SU(3)j is described by the Nambu-Goldstone mode, 7t°-{x) (a = 1, ■ • ■ , 8), 
and i2wz denotes the Wess-Zumino term. The ansatz for the chiral field of SO (3) soliton on 5*^ 
is represented with two profile functions, and ip{fi), asH, B] 



U = exp 



^{A-fx(/i)-(^(A-f)2-l)^(/x)} 



Here, A = (Ai,A2,A3) = (A7,— A5,A2) is the generators of the S0(3) group, where \a {a = 
1,2, •••,8) are Gell-Mann matrices, and three polar coordinates on are /i, 6, and 0. The 
baryon current density B^lfi) and the energy of the Skyrmion E{U) on with radius R can 
be easily obtained from those on -R'^H] by simple replacement! 



r ^ i? sin /i, 

dr ^ Rdfi, (3) 
dV = Anr'^dr dV = AnR^ sin^ /id/i. 

The winding number, which is identified with the baryon number, is 

B = AttR^ r sin^ /id/i5° (/i) , (4) 
Jo 
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where B^{fi) is the baryon number density, 



247r2 



■TilLaLjsL^ 



27r2i?3 sin^ /i 



dx.. 3 SdV' . .3 

-(l-cosxcos-V^) + --sinxsm-^ 



(5) 



with La = U'^daU ■ For the B = 2 SkjTrmion on 5''^, the profile functions satisfy the boundary 
conditions 

X(0) = 7r, tp{0) = ^n, x(vr)=0, V^(7r) = 0. (6) 
The static energy of the SO (3) Skyrmion on is written as 

EiU) = JdV &TTid^U^d''U)-^TT[d^UU\dMU^]' 



^ flRh + 



where Ii and I2 are dimensionless integrals, 
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n - cos X cos -ip 
sm /i I / 



(7) 



(8) 



2tt / d/i 
Jo 



— 2 — \ 3 sin^ X sin^ -ip + ( 1 — cos x cos -ip 
sin /i I 2 V 2 



+ ( 1 — cos X cos -ip 



'dx\\9_ /df\ ' 1 ^.^ ^ ^.^ 3 ^ dx # 
^d/i/ 4 yd/iy 1 2 d/x d/x 



(9) 



The Euler- Lagrange equations for x and ip are obtained by extremizing the static energy with 
respect to them. 



3 Numerical results and concluding remarks 

The profile functions are obtained numerically by extremizing the energy functional E{U) 
under the boundary conditions (H). We show in Fig.l the profile functions of the B = 2 SO (3) 
Skyrmion on in typical two cases using the parameters in ref. 0: /^r = 67 MeV and e = 4.47. 
(In Refs. [1T2|, |13|, the energy scale is rescaled by /^/(-\/2e) = 10.6 MeV, and the length scale 
by l/(-\/2e/^) = 0.466 fm to reduce variable dimensionless.) There are always two physically 
equivalent solutions: one is localized around the north pole, and another is localized around 
the south pole when R is large. The two solutions, (xi? "^i) and {x2, '^2), have the same energy, 
and have the relation, 



X2(/i) = 7r-xi(7r-/i), 
2 

^2(/^) = -71 -tpi{7f - fi). 



(10) 



At low densities, the Skyrmion localizes at either pole. In Fig. 1(a), the profile functions for 
i? = 2 fm are shown, where the Skyrmion is localized around the south pole. The behavior 
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of the profile functions is similar to that for the Skyrmion on the fiat space -R'^IQ. On the 
other hand, as shown in Fig. 1(b), the profile functions spread all over the hypersphere at high 
densities. 

Fig. 2 shows the energy of the SO (3) Skyrmion solution with the lowest energy as a function 
of the hypersphere radius R. At large R, the energy of the Skyrmion on approaches the one 
on R^. The behavior of the energy is essentially parallel with the SU(2) Skyrmion on S'^[0. 
The character of the lowest energy solution changes around i? ~ 1 fm. The minimum energy is 
attained at i? ~ 0.57 fm, and the minimum value is only 4 % above the BPS saturation value 
Ebps = 6v^7r^(/7r/e)|-B| (Bogomol'ny bound). It can be shown that the BPS saturation [|15 



IS 



never attained for the SO (3) Skyrmion on unlike the SU(2) Skyrmion on S^. In the B = 1 
SU(2) case, the BPS saturation is achieved by the "identity map" , f/(/i, 9, ip) = cos /(/x) + ir ■ 
fsin/(/x) with /(/i) = /i or /(/x) = tt — /i [|l^]. The SU(2) solution on 5*^ becomes "identity 



map" for small as R < Re- On the other hand, in the S0(3) soliton solution, there are 
two kinds of profile functions, x(/^) and ip{fi), and the B = 2 S0(3) solution on 5*^ tends to 
approach x(/i) = |'?/'(/i) = vr — /i as i? decreases as shown in Figs. 1(a) and (b). Moreover, there 
is a similarity in that the profile function has linear dependence on /i, so that we call the B = 2 
S0(3) solution with xif^) = fV'l/^) = vr — yU as the "identity map" in this paper. Being different 
from the SU(2) case, the "identity map" does not satisfy the Euler-Lagrange equations so that 
it is slightly above the lowest energy curve at high densities. 

At the low density phase, the baryon density is localized around either the north or the south 
pole, while at the high density phase the SO (3) Skyrmion spreads over the hypersphere, and 
loses its identity as a soliton. This may correspond to a phase transition from the confinement 
phase to the deconfinement phase. In Fig. 3, we show the baryon number densities B^{fi) of 
the solutions of several hypersphere radii around this phase transition. We can clearly see the 
localized-delocalized transition. 

The phase transition occurs so smoothly that it is difficult to find the critical density strictly. 
If we put the critical radius Rc to be 0.8~1.0 fm, the corresponding critical density pc = 
2/(27r^i?|?) is 0.1~0.2 baryons/fm^, which is the same order of magnitude as the equilibrium 
nuclear matter density po- This is rather smaller value than the expected critical density which 
is the order of several times po- Too low critical density problem appears also in B = 1 hedgehog 
Skyrmion matter |jTl], |T2|, |13[. The critical density pc is 0.07 ~ 0.18 fm~^ in these models. One 
of the reason is probably the lack of the kinetic energy as stated in ref . [0] . Nevertheless, it is 
interesting to be successful in modehng the high density phase transition in QCD qualitatively 
by such simple models. 

The order parameter ctq, which is a measure of the chiral symmetry breaking, is introduced 
by the spacial average [p!6|, 

/ dV-TiU = - r sin^ pdp ( -e'"^ + - cos xe"^ 
Js-^ 3 n Jo V3 3 



fill 



As shown in Fig. 4, the order parameter cxo becomes smaller at high densities, which is considered 
to be the sign of the chiral symmetry restoration. The non-zero value of the order parameter 
may be due to the limitation of the non-linear sigma model. 
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The Skyrmion root mean square (r.m.s.) radius is given by 



f = y^((i?sin/i)2) = ^ J^^ dy50(/i)(i?sin/i)2 ^ 2 (12) 

and shown in Fig.5. The contributions from the upper half (0 < /x < 7r/2) and the lower half 
(7r/2 < < 7r) of the hypersphere are also shown. At low densities, f is almost determined by 
the contribution from the half of the hypersphere where the Skyrmion localizes, while the both 
contributions are almost the same at high densities. Thus the deconfinement phase transition is 
read also in Fig.5. At low densities, the Skyrmion size coincides with the size of the Skyrmion on 
R^. In the localized phase {R> Rc — 0.8 — 1.0 fm), the size f of the S0(3) Skyrmion becomes 
larger corresponding to the spread of its profile functions for small R or at high densities. This 
enlargement implies the swelling of the H-dibaryon in high density H-dibaryon matter. In the 
deconfinement phase, where the Skyrmion loses its identity as the soliton, f decreases linearly 
with the hypersphere radius R. 

Finally, we compare the S0(3) soliton solution with the SO(2)xSO(2) one, which is also 
a candidate of the ground state of the B = 2 Skyrmion on a small hypersphere in the 
two-flavor case [jl3|. We show in Fig. 6 the total energy as a function of the hypersphere radius 
R for the S0(3) and the SO(2)xSO(2) soliton solutions in the B = 2 sector on S^. The energy 
difference is found to be quite small between the S0(3) and the SO(2)xSO(2) solutions for the 
high density region as i? ^ 1 fm. For instance, one finds 0.90 < -Eso(3)/-£'SO(2)xso(2) < 1-03 for 
0.2 < i? < 1.0 fm. One also finds Eso(3) > 1.041Ebps and Sso(2)xso(2) > 1.045Ebps. From the 
detailed analyses, the S0(3) solution is the true ground state for R > 0.99 fm (dilute case) and 
R < 0.64 fm (dense case), while the SO(2)xSO(2) solution is more stable for the intermediate 
density regions as 0.64 < R < 0.99 fm. Thus, the S0(3) soliton solution is the true ground 
state for the small hypersphere with R < 0.64 fm in the A^f = 3 chiral limit in the mathematical 
sense. In the real world, however, the large symmetry breaking in the strangeness sector should 
enhance -Eso{3) lai'gely, and therefore the SO(2)xSO(2) Skyrmion may become more stable. 

We next consider the chiral order parameter ao in eq. (pUD . Since the SO(2)xSO(2) solution 

takes the form as f/ = ( ^ j with j e'^^ '^^dV = |jl3|] independent of i?, ao = |. It is 

to be noted that the SO(2)xSO(2) solution on 5''^ provide the "chiral symmetric phase" only 
in the two-fiavor sector, and does not lead to the Af = 3 chiral symmetric phase. 

In summary, the SO (3) Skyrmion solution with B = 2 has been investigated on a hy- 
persphere 5''^ to study high density H-dibaryon matter. In this system, the deconfinement 
transition and the chiral symmetry restoration have been found at high densities similarly to 
B = 1 SU(2) Skyrmion matter. As a remarkable fact, the H-dibaryon in H- matter swells at 
high densities. The swelling phenomenon is considered to refiect the attractive nature of the 
interaction between H-dibaryons[|l7l. Such a possibility of the swelling of the H-dibaryon would 
be important for the study of H-matter in dense neutron stars, because the compactness of the 
H-dibaryon is a basic assumption for the existence of H-matter. For the realistic treatment of 
H-dibaryon, we need to take account of the SU(3) flavor symmetry breaking. This will be done 
in the future work. 
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Figure 1: The profile functions, and ip{fi), of the B = 2 S0(3) Skyrmion on vs. the 
angle /i of the hypersphere for (a) R — 2.0 fm, and (b) R — 0.1 fm. 
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Figure 4: The change of the chiral order parameter ctq with respect to the hypersphere radius 
R. 
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Figure 5: The Skyrmion r.m.s. radius f defined by eq.(|T^) as a function of tlie liyperspliere 
radius R (solid line). The contributions from the upper half (0 < fi < 7r/2) and the lower 
half (7r/2 < /i < tt) of the hypersphere are also shown by dashed hne and dot-dashed hne, 
respectively. 
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Figure 6: The comparison between the S0(3) and the SO(2)xSO(2) solutions on . The 
energies, ii^so(3) and -Eso(2)xso(2), are plotted as functions of the radius R. 
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